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Amendments of the Claims: 

A detailed listing of all claims in the application is presented below. This listing of claims 
will replace all prior versions, and listings, of claims in the application. All claims being 
currently amended are submitted with markings to indicate the changes that have been made 
relative to immediate prior version of the claims. The changes in any amended claim are being 
shown by strikethrough (for deleted matter) or underlined (for added matter). 

Claims 1-14. (Canceled) 

15. (Currently Amended) A practical numerical computer-implemented method for reliably 

computing a dynamical decomposition point of a stable equilibrium point for large-scale 
nonlinear systems, comprising the steps of: 

a) given a stable equilibrium point x s; 

b) moving along a search path p t (x s ) = {x s + * x *, t e 91 + } starting from x s and 

detecting an exit point, x^, at which said search path (p t (x s ) exits a 
stability boundary of a stable equilibrium point x s ; 

c) using said exit point x„ as an initial condition and integrating a nonlinear 

system (4t3> to an equilibrium point x d ; and 

d) computing said dynamical decomposition point with respect to the stable 

equilibrium point x^ wherein said search path is x d .; and 

e) displaying the dynamical decomposition point. 

16. (Currently Amended) The method of claim 15, wherein a method for computing said exit 

point of the nonlinear system (4r3) comprises the step of moving along said search path 
9t a {*., + * x S, t G 91 + } starting from x s and detecting said exit point x a , which is 
a first local maximum of an objective function C(x) along said search path <p, (x s ) . 
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17. (Currently Amended) The method of claim 15, wherein a method for computing a dynamical 

decomposition point comprises the steps of; 

a) using said exit point x €X as an initial condition and integrating a nonlinear 

system fkO) to a first local minimum of a norm ||F(x)|j along the 
corresponding trajectory, where F{x) is a vector field of the nonlinear 
system 1 4=3}. and letting the point at which the first local minimum of 

occurs be denoted x° d , and is called the minimum distance point 
(MDP); and 

b) using said MDP x° as an initial guess and solving a set of nonlinear algebraic 

equations of said vector field (4r3> F(x) = 0, wherein a solution is x d , and 
a dynamical decomposition point with respect to the local optimal solution 
x a and said search path <p t (x s ) is x d , 

18. (Currently Amended) The method of claim 15, wherein a method for computing said exit 

point with respect to a stable equilibrium point of the nonlinear system {4t2) and a search 
vector comprises the step of computing an inner-product of said search vector and the 
vector field of the nonlinear system £k£) at each time step, by moving along said search 
path <p t (x s )={x s + txs, t e 91 + } starting from x s and at each time-step, computing an 
inner-product of said search vector $ and vector field F(x) y such that when a sign of 
said inner-product changes from positive to negative, said exit point is detected. 

19. (Currently Amended) The method of claim 15, wherein a method for computing said exit 

point of thenonlinear system <43) with respect to a stable equilibrium point and a search 
vector comprises the steps of: 

a) moving from said stable equilibrium point along said search vector until an 
inner-product of said search vector and the vector field of the nonlinear 
system {4£)-changes sign between an interval [f r , t 2 ]\ 
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b) applying a linear interpolation to an interval [t u r 2 ]> which produces an 

intermediate time t 0 where an interpolated inner-product is expected to be 
zero; 

c) computing an exact inner-product at fa such that if said value is smaller than a 

threshold value, said exit point is obtained; and 

d) if said inner-product is positive, then replacing t\ with fa and otherwise 

replacing t 2 with / 0 and going to step b). 

20. (Currently Amended) The method of claim 15, wherein a method for computing a minimum 

distance point (MDP) of the nonlinear system (4.2) satisfying required c onditions4Cj4 
and (C2) comprises the steps of: 

a) using said exit point as an initial condition and integrating the nonlinear system 

(4t3) for a few time-steps, and letting the end point be denoted as the 
current exit point; 

b) checking convergence criterion, and, if a norm of said current exit point 

obtained in step a) is smaller than a threshold value, then declaring said 
point as said MDP and stopping the process, otherwise, going to step c); 
and 

c) drawing a ray connecting a current exit point on a trajectory and a local optimal 

solution (equivalently, a stable equilibrium point), replacing said current 
exit point with a corrected exit point, which is a first local maximal point 
of objective function along said ray, starting from the stable equilibrium 
point, and assigning this point to said exit point and going to step a). 

21 . (Currently Amended) The method of claim 15, wherein a method for computing said 

dynamical decomposition point of the nonlinear system (4t2) satisfying required 
conditions (CI) and (C2) with respect to a stable equilibrium point x s and a search vector 
$ , comprises the steps of: 
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a) moving along said search path <p t (x x ) s {x s + t x s,t £ 91 + } starting from x x and 

detecting a moment that an inner-product of said search vector s and the 
vector field F(x) of thejionlinear systemn(4T3) changes sign, between an 

interval [f p f 2 ], stopping this step if t t is greater than a threshold value and 
reporting that there is no adjacent local optimal solution along this search 
path, otherwise, going to step b); 

b) applying linear interpolation to said interval [t v t 2 ] > which produces an 

intermediate time / 0 where said interpolated inner-product is expected to 

be zero, computing an exact inner-product at t 0 , and if said value is 

smaller than a threshold value, said exit point is obtained, and going to 
step d); 

c) if said inner-product is positive, then replacing with t 0 , and otherwise 

replacing t 2 with t 0 and going to step b); 

d) using said exit point as an initial condition and integrating a the nonlinear 

system (4^3) for a few time-steps, and letting the end point be denoted as 
the current exit point; 

e) checking convergence criterion, and if a norm of said point obtained in step d) 

is smaller than a threshold value, then declaring said point as the MDP and 
going to step g), otherwise going to step e); 

f) drawing a ray connecting a current exit point on said trajectory and a stable 

equilibrium point, replacing said current exit point with a corrected exit 
point which is a first local maximal point of objective function along said 
ray starting from the stable equilibrium point, and assigning this point to 
said exit point and going to step d); and 

g) using said MDP as an initial guess and solving a set of nonlinear algebraic 

equations of the vector field of the nonlinear system {4r3) F(x) = 0 , 
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wherein a solution is t d , such that said DDP with respect to a stable 
equilibrium point x s and search vector s is x d . 

22. (Cancelled) 

23. (Currently Amended) A hybrid search computer-implemented method for obtaining a local 

optimal solution of a^eneral unconstrained nonlinear programming problem {4r+} 
starting from any initial point, comprising the steps of: 

a) given an initial point jc 0 ; 

b) integrating a nonlinear dynamical system d e scrib e d by (4.2) that satisfies 

required conditions (CI) and (G2 - ) from said initial point x 0 to obtain a 

trajectory fa [x 0 ) for n time-steps, n being an integer, and recording the 
last time-step point as the end point, and if it converges to a local optimal 
solution, then stopping, otherwise, going to step (c); 

c) monitoring a desired convergence performance criterion along the trajectory 

fa (jto) in terms of the rate of decreasing values in the objective function 
under study, and if the desired criterion is satisfied, then using the end 
point of trajectory fa(x 0 ) as the initial point and going to step (b) 7 
otherwise, going to step (d); and 

d) applying an effective local optimizer (i.e., a method to find a local optimal 

solution) from said end point in step (b) to continue the search process, 
and if it finds a local optimal solution, then displaying the local optimal 

solution and stopping, otherwise, setting the end point of trajectory fa(x 0 ) 

as the initial point, namely x 0 , and going to step (b). 

24. (Cancelled) 
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25. (Previously Presented) The method of claim 44-23, wherein in step (a) the initial point * 0 is 

given by a method comprising the steps of: 

a) moving along a search path starting from a local optimal solution x opt and 

applying a DDP search method to compute a corresponding DDP, and if a 
DDP can be found, then going to step (b), otherwise, trying another search 
path and repeating this step; 

b) letting said DDP be denoted as x d , and if x d has previously been found, then 

going to step a), otherwise going to step c); and 

c) setting x Q = x^ + (1+ e)(x d - x opt ) where £ is a small number, and applying a 

hybrid search method starting from ;e 0 to find a corresponding adjacent 
local optimal solution. 

26. (Previously Presented) A computer-implemented method for obtaining a global optimal 

solution of unconstrained nonlinear optimization problems, comprising the steps of: 

a) choosing a starting point; 

b) applying the hybrid search method of claim 23 using said starting point to find 

a local optimal solution x* ; 

c) setting; = 0, V s = , Vj ew = and V d = {<* } ; 

d) wherein setV^ 1 = {<f>} and for each local optimal solution in (i.e., */), 

performing steps (e) through (k); 

e) defining a set of search vectors Sf , i = 1 ,2, . . . , m h and setting i = 1 ; 

f) if i > rrij , then going to step (1); otherwise, applying a DDP search method 

along the search vector 5/ to find a corresponding dynamical 
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decomposition point (DDP), and if a DDP is found, then going to step (g), 
otherwise, setting i = i + 1 and going to step (f); 

g) letting the found DDP be denoted as x J di , checking whether it belongs to the 

set V d , i.e., x J di e V d ?, and if it does, then setting i = i + 1 and going to step 
(f), otherwise settinjg V d = V d v {jcj jand going to step (h); 

h) for the DDP x J di , performing steps (i) through (j) to find a corresponding 

adjacent local optimal solution; 

i) letting - x J s + (l + - */ ) , where e is a small number; 

j) applying said hybrid search method using as the initial condition to find the 
corresponding local optimal solution, and letting it be denoted as x i sj ; 

k) checking whether y has been found before, Le., x^ e V x ?, and if it has 

already been found, then setting i = / + 1 and going to step (f), otherwise, 
setting V s = V s u {<;} and v£ = V^ 1 u and setting f = i + 1 and . 
going to step (f); 

1) examining the set of all newly computed local optimal solutions, V*£ » and if 
Vj£ is non-empty, then setting j= j + I and proceeding to step (d), 
otherwise proceeding to the next step; and 

m) identifying the global optimal solution from said set of local optimal solutions 
V s by comparing their corresponding objective function values in set V, 
and; 

n) displaying the global optimal solution. 
27. (Cancelled) 
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28. (Currently Amended) A computer-implemented method for obtaining a global optimal 

solution of a constrained nonlinear programming problem (4r§), comprising the steps of: 

a) Phase I: finding all feasible components of the constrained nonlinear 

programming problem wherein one effective local method is combined 
with said dynamical trajectory method, comprising the following steps to 
find a feasible solution of the constrained optimization problem (4r7): 



i) given an initial point *o; 



ii) integrating a nonlinear dynamical system d e scrib e d by (4.9) that 

satisfies required c onditions (CI 1) and (CI 2) from said initial 




and recording the last time-step point as the end point, and if it 
converges to a feasible solution, then stopping, otherwise, going to 
step (iii); 



iii) monitoring a desired convergence performance criterion in the 

objective function IM^i, a vector norm of H ^ in^4^ the 
constrained optimization problem , and if the desired criterion is 




point and going to step (ii), otherwise, going to step (iv); and 



iv) applying an effective method from the said end point in step (b) to 
continue the search process, and if it finds a feasible solution of 
(4^£) the constrained optimization problem then stopping, 




the initial point and going to step (ii) ; and 



b) Phase II: finding all local optimal solutions for the constrained nonlinear 
programming problem in each feasible component: 
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c) choosing a global optimal solution for the constrained nonlinear programming 

problem from the local optimal solutions found in step (b); and 

d) displaying the global optimal solution for the constrained nonlinear 

programming problem. 

29. (Currently Amended) The method of claim 28, comprising the steps of: 

a) approaching a path-connected feasible component of a constrained 

optimization problem-£k5); and 

b) escaping from said path-connected feasible component and approaching 

another path-connected feasible component of said constrained 
optimization problem-(4r§). 

30. (Currently Amended) The method of claim 28, comprising the steps of: 

a) approaching a stable equilibrium manifold of a nonlinear dynamical system 

(4r9) satisfying required c onditions- ft51 - n and (C\ - 2) ; and 

b) escaping from said stable equilibrium manifold and approaching another stable 

equilibrium manifold of said nonlinear dynamical system (4r9) satisfying 
said condition s (CI 1) and (CI - 2) . 

31. (Currently Amended) The method of claim 28, comprising the steps of: 

a) starting from a point in a feasible component and approaching a local optimal 

solution located in said feasible component of an optimization problem 
<4^?); and 

b) escaping from said local optimal solution and approaching another local 

optimal solution of said feasible component of said optimization problem 

32. (Currently Amended) The method of claim 28, comprising the steps of: 
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a) in a deterministic manner, first finding all stable equilibrium manifolds of a 

nonlinear dynamical system that satisfies required c onditions (CI 1) and 
(Gi-2); 

b) in a deterministic manner, finding all stable equilibrium points of a nonlinear 

dynamical system that satisfies sakLconditions (C2 1) and (C2 2) ; and 

c) then from said stable equilibrium point, finding a global optimal solution. 

33. (Currently Amended) The method of claim 32, comprising the steps of: 

a) given a feasible solution of constrained nonlinear programming problem-(4T#); 

b) finding a stable equilibrium point of a nonlinear dynamical system frl -r lO) 

satisfying required c onditions fC2 1) and (C2 2) ; 

c) moving from said stable equilibrium point to a dynamical decomposition point, 

in order to escape from a local optimal solution; and 

d) approaching another stable equilibrium point of said nonlinear dynamical 

system (4r l -Q- ) satisfying said c onditions (C2 1) and (C2 2) in the same 
path-connected feasible component, via said dynamical decomposition 
point. 

34. -36. (Cancelled) 

37. (Currently Amended) A hybrid search method for Phase II of claim 28, comprising the 

following steps to find a local optimal solution of the constrained optimization problem: 

a) given a feasible point x 0 ; 

b) integrating a nonlinear dynamical system described by (1.10) that satisfies 

required c onditions ( G2 1) and (C2 - 2 )-from said initial point to obtain a 
trajectory tf,(x 0 ) for n time-steps, n is an integer, and recording the last 
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time-step point as the end point, and if it converges to a local optimal 
solution, then stopping, otherwise, going to step (c); 

c) monitoring a desired convergence performance criterion in an objective 

function C(x ) of (4.5). , and if the desired criterion is satisfied, then using 
the end point of trajectory $ t ) as the initial point and going to step (b), 
otherwise, going to step (d); and 

d) applying an effective method from the said end point in step (b) to continue the 

search process, and if it finds a feasible solution of (434 the constrained 
optimizati on problem then stopping, otherwise, setting the end point of 
trajectory $ t (x 0 ) in step (b) as the initial point and going to step (b). 

38. (Currently Amended) The method of claim 30, comprising the steps of: 

a) choose a starting point; 

b) initialization j = 0; 

c) applying a hybrid search method for Phase I using said starting point to find a 

feasible point in a (path-connected) feasible component, setting a point so 
found as an initial point, and applying the hybrid search method for Phase 
II to find a local optimal solution x/ 0 ; 

d) starting from said initial point x J s 0 , applying a numerical method for Phase II to 

find all local optimal solutions in said feasible component and recording 
them as the set VJ and set V s ~ V s U Vj ; 

e) for the local optimal solution x/ 0 , defining a set of search vectors 

5/,/ =1,2,...,*,.; 

f) if i > then going to step (1), otherwise, going to step (g); 
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g) for each search vector Sf , applying a reverse-time trajectory method to find a 

point lying in an unstable equilibrium manifold of a system {4t9)> and if a 
point is found, letting it be denoted as x' dJ and going to step (h), otherwise 
setting / = i + 1 and going to step (f); 

h) setting x' 0j = x J s + 0* + £)(x l dj - x J J , where e is a small number and x J s is a local 

optimal solution selected in step (e), and applying the hybrid search 
method for Phase II using x' 0 J as the initial point to find a point lying in a 
stable equilibrium manifold of system (4r9) satisfying required conditions 
(CI 1) and (CI 2) , and letting the solution be denoted as x° . : 

S,J 

i) starting from said initial point x* %J and applying the hybrid search method of 

step (b) of claim 28 to find a local optimal solution in said (path- 
connected) feasible component jc^ 1 ; 

j) checking whether x/* 1 has been found before, Le. x J s + l e V x ?, and if it has been 

bound before (Le., the said feasible component has been visited before), 
then setting i - i + 1 and going to step (f), otherwise setting 

V. = K u {*io} and V£ = v£ u {xtf } and going to step (k); 

k) starting from said initial point x^ , applying a numerical method for Phase II 
to find all local optimal solutions in said feasible component and recording 
them in the set Vj +l , set V 5 = V, u V/ +1 and setting i = i -h 1 and going to 
step(f); 

1) examining the set of all newly computed local optimal solutions V^J , and if 
V£ v l is empty, then going to the next step, otherwise, setting; =7 + 1 and 
going to step (e); and 
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m) identifying the global optimal solution from said set of local optimal solutions 
V s by comparing their objective function values. 

39. (Currently Amended) A numerical method for Phase II of claim 28, comprising the following 
steps to find all the local optimal solutions in said path-connected feasible component: 

a) given an initial point lying in said path -connected feasible component; 

b) initialization; 

c) applying a hybrid search method for Phase II using the initial point to find a 




d) for each local optimal solution in Vj^ (i.e., xj ), performing steps (e) through 



00; 



e) defining a set of search vectors Sf J = l,2,...,/n ; . , and setting i = 1; 



f) if i > m h then going to step 0), otherwise, applying a DDP search method, 



wherein a nonlinear dynamical system £MO)-satisfies required conditions 
(C2 1) and (C2 2) to find a corresponding dynamical decomposition point 
(DDP), and if a DDP is found, then going to step (g), otherwise, setting i - 
i + 1 and going to step (f); 



g) letting the found DDP be denoted as xj . , checking whether it belongs to the 

set V d , i.e., xj , e V d ?, and if it does, then setting i = i + 1 and going to step 




h) for the DDP jcj ; , performing steps (i) through (j) to find a corresponding 
adjacent local optimal solution; 
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i) setting x^ = xf + (l + £){x^ - ) , where 6 is a small number; 

j) applying a hybrid search method using *£. as the initial condition to find the 
corresponding local optimal solution, and letting it be denoted as xf • ; 

k) checking whether x f sj has been found before, Le. x l sj e V s ?, and if it has 

already been found, then setting i= i + 1 and going to step (f), otherwise, 
setting V, = V s u {<,] and v£ l = V£ u {<,} and setting i = £ + 1 and 
going to step (f); and 

1) examining the set of all newly computed local optimal solutions Vj£ n , and if 
v £Zm is non-empty, then setting,/ = j +1 and proceeding to step (d), 
otherwise outputting all the local optimal solutions in said path-connected 
, feasible component contained in set V s and stopping the process. 

40. (Previously Presented) A computer-implemented method for obtaining the global optimal 
solution of constrained nonlinear programming problems, comprising the steps of: 

a) using a transformation technique to transform a constrained optimization 

problem into an unconstrained optimization problem, then applying the 
following steps to find the global optimal solution of the unconstrained 
optimization problem; 

b) choosing a starting point; 

c) apply a hybrid search method using said starting point to find a local optimal 

solution jc°; ■ 

d) setting; = 0, V, = , = and V H = {^} ; 

e) setting Vj£ = and for each local optimal solution in (i.e., 

performing steps (e) through (k); 
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f) defining a set of search vectors S? , i = 1,2, . . m h and setting i = 1; 

g) if f > m jy then going to step (1), otherwise, applying the DDP search method of 

claim 15 along the search vector S? to find a corresponding dynamical 
decomposition point (DDP), and if a DDP is found, then going to step (g), 
otherwise, setting i = / + 1 and going to step (f); 

h) letting the found DDP be denoted as x J d . , checking whether it belongs to the 

set V d , le.i x J d . e V d ?, and if it does, then setting x = i + 1 andgoing to step 
(0, otherwise setting V d = V d u (xj ( .} and going to step (h); 

i) for the DDP x J dJ , performing steps (i) through (j) to find a corresponding 

adjacent local optimal solution; 

j) setting = xj + (l+ e)[x J d . - */) , where e is a small number; 

k) applying a hybrid search method using as the initial condition to find the 
corresponding local optimal solution, and letting it be denoted as x l s } ; 

I) checking whether x' sj has been found before, i.e. x l s j e V s ?, and if it has 

already been found, then setting / = i + 1 and going to step (f), otherwise, 
setting V s = V x u {<,} and v£ ! = u {< .} and setting i = r + 1 and 
going to step (f); 

m) examining the set of all newly computed local optimal solutions, Vj£ , and if 
v new is non-empty v then setting j= j+l and proceeding to step (d), 
otherwise proceeding to the next step; aad 
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n) identifying the global optimal solution from said set of local optimal solutions 
V s by comparing their corresponding objective function values in set V s ; 
and 

o) displaying the global optimal solution. 
41. (Cancelled) 
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